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Abstract 

The decay B — > J/ipKs remains the most important source of information for the Bd 
mixing phase, determined by the CKM angle /3 in the standard model. When aiming at 
a precision appropriate for present and coming high luminosity colliders, the corresponding 
hadronic matrix elements are a major obstacle, as their precise calculation is still not fea- 
sible with existing methods. Flavour symmetries offer a possibility to extract them from 
data, however again with limited precision. In this article, we propose a framework to take 
subleading contributions in B lh d, s — ► J/tpP decays into account, P £ {ir,K, {rjs)}, using an 
SU(3) analysis, together with the leading corrections to the symmetry limit. This allows for 
a model-independent extraction of the Bd mixing phase adequate for coming high precision 
data, and additionally yields information on possible New Physics contributions in these 
modes. We find the penguin-induced correction to be small, |A5| < 0.01, a limit which can 
be improved with coming data on CP asymmetries and branching ratios. Finally, the sensi- 
tivity on the CKM angle 7 from these modes is critically examined, yielding a less optimistic 
picture than previously envisaged. 
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1 Introduction 



With the LHC up and running, we have entered a new era in the search for new physics (NP), 
as well as in the determination of fundamental parameters of the standard model (SM). In 
flavour physics, a central role is played by the weak phase entering the CKM matrix, and the B 
factories together with the Tevatron experiments have done an amazing job in confirming the 
CKM mechanism as the main source of low-energy CP violation. The 0(%)-measurement of 
the CKM angle j3 has been a highlight in that respect. This precision became possible due to 
the fact that in the "golden mode", Bj — > J/ipKs, explicit calculation of the relevant matrix 
elements can be avoided, once subleading doubly Cabibbo suppressed terms are assumed to 
vanish jl|, in combination with a final state with a very clear signature. However, given the 
precision the LHC experiments and planned next-generation B factories are aiming at for this 
mode and related ones, a critical reconsideration of the assumptions used is mandatory. 

The endeavour of including subleading contributions has been started already some time 
ago, see e.g. [2}j9]. As pointed out there, in addition to including these terms, the procedure 
allows for an improved access to NP in B meson mixing, as the SM correction constitutes a "fake 
signal" in that analysis. This is especially interesting, as for several years different extractions 
of the corresponding angle /3 in the unitarity triangle (UT) have been only marginally consistent 



(see e.g. 10-14]), especially when including B — > tv in the fit. 

The main problem lies in the evaluation of hadronic matrix elements, which still does not 
seem feasible to an acceptable precision for the decays in question. To avoid explicit calculation, 
typically symmetry relations are used, which again have only limited precision. Furthermore, 
most of the existing proposals to address this question use (mild) additional dynamical assump- 
tions. In this work, we address both issues: we use symmetry under flavour SU(3) to relate 
B —7- J/vpP decays, where P stands for a pseudoscalar (octet) meson, without neglecting the 
subleading contributions; in addition, we model-independently include corrections to this limit, 
without referring to factorization assumptions. Doing so, we demonstrate the importance of 
SU (3)-breaking contributions in this procedure, whose neglect can lead to an overestimation 
of the subleading effects. Furthermore, the additional assumption used in [3j|5] is discussed in 
a more quantitative manner, to be able to determine a possible influence from this source as 
well. Apart from these points, the decays in question also have some sensitivity to the CKM 
angle 7, see e.g. [2j[8]. While these modes are not competitive in that respect to e.g. B — > DK, 
the extraction could serve as an independent crosscheck, and in the event of a largely deviating 
value even as a hint of NP in the decay amplitudes. 

The outline of this article is as follows: we introduce B — > J/ipP decays in section [2j 
pointing out the main features of the decay amplitudes which facilitate the following analysis, 
and presenting the available data. In section |3j we perform the SU(3) analysis, including linear 
breaking terms. The resulting expressions are then used in section [4] to perform first fits to the 
available data, to illustrate the impact future measurements are going to have, and to investigate 
the potential influence of NP. We conclude in section [5j 



2 Non-leptonic B — > J/i/?P decays 

The main focus of this work lies on the decay Bj — > J/ipK, from whose indirect CP asymmetry 
the Bd mixing phase can be extracted. It is related by flavour symmetry to the group of modes 
B — > J/tpP, where B E {B~ , B° , B s }, and P denotes a pseudoscalar meson from the same 
multiplet as the kaon, i.e. P G {K , n , (rjs)} . Using the unitarity of the CKM matrix, the 
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relevant SM hamiltonian has the structure (see e.g. |15| and references therein) 

n 



p=d,sq=u,c \i=l i=3 / 



= Y, n ^ p + n ^ p ^ n c + n u , (2.1) 

p=d,s 

where Oi are four-quark operators stemming from tree (i = 1,2), penguin (i = 3 . . . 6) and 
electroweak penguin (i = 7 . . . 10) diagrams. It describes both classes of decays relevant here, 
corresponding to a change in strangeness AS = — 1 or AS = 0. The corresponding b — > s tran- 
sitions are famous for being theoretically clean [lj, i.e. to very good approximation dominated 
by one hadronic amplitude which cancels in the expression for the indirect CP asymmetry. Es- 



timates yield corrections of the order 0(1O -3 ), only 16-18 ; it is however notoriously difficult 
to actually calculate the relevant matrix elements, and non-perturbative enhancements cannot 
be excluded. Given this fact, and the impressive experimental precision the flavour factories 
and the Tevatron have reached already, and the LHC and future Super Flavour Factories (SFF) 
are about to achieve, subleading contributions should be taken into account [2]-[9]. To this aim, 
the above mentioned symmetry-related processes are important, especially those with b — >• d 
transitions. For them, the suppression of the subleading contributions is much weaker, as the 
relative CKM factor is lacking the A 2 suppression, | V^Kdl/l Kb^edl ~ R u = y/p 2 + fj 2 ~ 0.37. 
This gives a handle to constrain them by data directly, making the extraction of the B& mixing 
phase (j) possible without dynamical assumptions. However, the main amplitudes in these tran- 
sitions are Cabibbo-suppressed compared to their b — > s counterparts, making the experimental 
information scarce. Nevertheless, the high statistics measurements becoming available make a 
future precision analysis feasible. 

To fully exploit the available information, the approach proposed here is to make a full SU (3) 
analysis, and thus to include 5-6 different decay mode^] This has the additional advantage to 
allow for the model-independent inclusion of the leading SU (3)-breaking contributions, albeit 
at the price of a slightly more complicated analysis. These contributions turn out to be crucial 
for a reliable extraction of the correction we are aiming at. Finally, the extraction of the CKM 
angle 7 proposed in [2] is included and generalized automatically, such that we can investigate 
its precision in the presence of S'C/(3)-breaking corrections. While different aspects of this 



anaysis have been discussed to more or less extent in the past [2j|5j[7J[8 20 , the combination and 
extension proposed here is critical to achieve the desired precision for the decays in question in 
a model-independent manner. 

The experimental data we are going to use are given in Tab. [I] Specifically two recent 
measurements are noteworthy: The ratio of branching fractions 

_ BR(B S -> J/^K S ) 
RKK = BR(B° -+ J/^Ks) (2 - 2) 



has been measured by the CDF- and LHCb-collaborations [24 [25 ; their values agree well with 
each other, leading to the average given below. This provides the first input from a B s decay to 
this analysis. Furthermore, the ratio 

_ BR(B- -> JjjjTT) 
R * K = BR(B- ^ J/jK=) (2 - 3) 



1 We will consider modes including singlets, i.e. Bj, s — > J/iprf in a future publication 19 
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Decay 


BR/W~ A A CP /% S CP Ref. 


5- -»• J/V'vr- (WA) 
dataset 2 (LHCb) 
5 s -»■ J/V'K 


8.71 ±0.32 0.6 ± 2.1 0.665 ± 0.022 |2l[ 
0.176 ±0.016 10 ±13 -0.93 ±0.15* 21 
10.13 ±0.34 0.1 ±0.7 2: 
0.50 ±0.04 1±7 2: 
0.39 ±0.02 0.5 ±2.9 % 
0.33 ±0.04 (24 


22 

22 

L 

L 

\ 

25 



Table 1: Experimental data used for the analysis. The values used from HFAG are from after the 
Moriond 2012 update. For B~ — > J/ipTr~ the first line corresponds to the former world average data 
(dataset 1), the second to the new LHCb results (dataset 2), as explained in the text. * This value is 
taken from HFAG, without enhancing the error due to differences of the two averaged measurements. 



has been measured by the LHCb collaboration |23| about 3 standard deviations away from the 
former world average of all direct measurements of this quantity. Instead of averaging the two 
results we will perform fits with each of them individually to quantify their influence. We call 
dataset 1 the one including the previous world averages for R^k and Aqp(B~ — > J/ipir~), and 
dataset 2 the one with the recent measurements of LHCb instead. 

Apart from these recent developments, the most striking feature of the data is the compa- 
tibility of all direct CP asymmetries with zero, constraining mainly the imaginary parts of the 
introduced strong amplitudes. For the CP asymmetries, we use the notation 



a C p(t) 



r(B d>s -> J/^X°) - T(B d)S -> J/ifjX°) 

. _v Jlih X°\ p<W Am. j 

(2.4) 



T(B d>t> -> J/^X0) + T(B d)S -)> J/^X ) 

S(B d;S -» J/ipX Q ) sm(Am dtS t) + A C p(B d>a -» J/jiX ) cos(Am d;S t) 
cosh(Ar^t/2) - A ATda (B dj8 -> J/VX )sinh(Ar d)S V2) 



where A = {Ks,l^ }> S(B diS — > J/t/jX°) and Acp(B dtS — > J/i/jX ) are referred to as indirect 
and direct CP asymmetry, respectively, and S(B djS ->■ J/ipX ) 2 + Acp(B d)S — > J/tpX ) 2 + 
AAT d , a ( B d,s -> J/i>X ) 2 = 1 holds as well as AL d « 0. 

Another interesting observation is the relatively large central value of the rate difference 

„ K T(B° -»■ JM°) - f (5- -> J/ibK-) , . 

AfiE=^ V - n ( ^ Vt r ^ -0.037 ±0.025. (2.5) 

7 L(S° -> J/V'A" ) + r(B- -»• J/ipK-) V ; 

While this measurement is still compatible with zero, a confirmation of this central value would 
imply much larger contributions from isospin-changing (A/ = 1) operators than expected in the 
SM. 

A more extensive discussion of these observables and their relation to the parameters in 
question is given in Sec. |4| Note that, although the branching ratios are given explicitly in 
Tab. [TJ we use in the numerical analysis the measured ratios where available, in order not to 
correlate the data additionally. We use the values R KK = 0.038 ± 0.004, = 0.049 ± 0.004, 
and R s £k = 0.0383 ± 0.0013. For Rrk we apply the correction due to the different definitions 
typically used in theory and experiment, recently discussed in [26 and also applied in 25 . 



2 In the extraction of the correction, the authors use exact SU{3) symmetry. However, the assigned error is 
large enough to accommodate also deviations from that limit. 
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3 SU(3) analysis 



In this section, the expressions for the decay amplitudes including penguin and leading first 
order S'C/(3)-breaking contributions are given. We concentrate here on the main results of that 
analysis; technical details are relegated to appendix [A} 

We start by performing the analysis in the SU(3) limit. For % c , we observe that neglecting 
the AI = 1,3/2 part of the electroweak-penguin operators, the flavour structure amounts to a 
pure triplet, 

(3)o,o,-2/3, and 7^ d ~(3) 1/2i _ 1/2il/3 , (3.6) 

using the notation i z ,y for a representation R with isospin / and hypercharge Y. The 
neglected contributions are suppressed by small Wilson coefficients, and additionally by the 
fact that these matrix elements do not receive factorizable contributions, which we will call 
"dynamical suppression" in the following. We deem the combined suppression to be strong 
enough to neglect this contribution, even when considering S'C/(3)-breaking contributions. 

For H u , the situation is complicated by the AI = 1,3/2 part of the occuring tree operators, 
which have large Wilson coefficients and are only weakly CKM suppressed in b — )■ d decays. 
They are however again dynamically suppressed. With these terms, the SU(3) analysis yields 
three independent reduced matrix elements. We confirm the expressions obtained in 



20 . 


For 


1 16 


{18] 



correspond to a ratio r pcn < 8% of these matrix elements with the leading one. However, as we 
allow conservatively for their enhancement, we will consider ratios up to 50%, and in some cases 
even larger enhancements for illustration purposes. 

In B s — > J/ipTT°, the otherwise leading part of the amplitude vanishes due to isospin conserva- 
tion; the only contributions to this decay are from so-called penguin-annihilation and exchange 
topologies, which are expected to exhibit additional suppression ~ Aqcd/^6, and are usually 
neglected (see e.g. [3j|5]). We discuss the possible inclusion of these terms in Sec. [4} 

As the next step we include linear SU(3) breaking due to the finite light-quark mass dif- 
ferences. Neglecting isospin breaking, the relevant tensor structure is (8)o,o,o> driven by the 
strange-quark mass term. The tensor product of the leading-order hamiltonian T~L C with this 
breaking term leads to three more reduced matrix elements. These include implicitly the expan- 
sion parameter expected for this kind of analysis, e ~ m s /A ~ 20 — 30%, where A denotes some 
hadronic scale, as reflected e.g. in the decay constant ratios (fx/ fn — 1) ~ (/s s //e d — 1) ~ 20%. 
Without further dynamical information, this is also the estimate for their ratios with the leading 
matrix element, characterized by rgjj^ as defined in the appendix. To be conservative, we allow 
for rs{/( 3 ) ~ 40% in the following fits (60% for illustration). Note that r su ^ ~ 1 corresponds 
to the absence of any hierarchy, i.e. the reduction of SU(3) symmetry to isospin. There is no 
sign in B meson decays for such a drastic behaviour, and the fits discussed later indicate in fact 
a value of the expected order. For A u , we can neglect these corrections, as already the leading 
terms are strongly suppressed. To our knowledge, these decompositions are derived here for the 
first time, only a subset has been considered in (T 27 



As A C (B S — > J/^7T°) = was enforced by isospin at leading order, all breaking terms vanish 
as well, as we do not break isospin symmetry. The coefficients for charged and neutral B — > J/ipP 
decays (P £ {ir,K}) remain pairwise forced to be equal by isospin as well (modulo a factor of 
\/2 for the 7T°). Related to that is the fact, that the rank of the relevant coefficient matrix is 
not maximal: all amplitudes can be expressed by only two independent combinations, without 
any approximation, thereby reducing the number of linearly independent amplitudes needed to 
describe these decays from seven to six. 
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The framework developed so far yields a description of the six relevant decay amplitudes in 
terms of 11 hadronic parameters and the CKM phase 7, which reads 



A{B° -> J/4>K°) = M [1 + 2R €l + \ 2 e-^ (R ul + 5)] , 

A(B° -> J/^vr°) = -A^= [1 - i? el - i? e2 - e*i (Rui - 5)] , 

A(B~ -> J/rf>K-) = JV[l + 2R tl + X 2 e-^ (R u2 -5)] , 
A(B~ -> J/^vr") = -AAA [l - i? el - i? e2 - (R u2 - 5)] , 
A{B S -> J/^K°) = -AAA [1 - R el + R e2 - (R ul + 5)} , 

A(B S ^ J/4>tt°) = ^-~X 2 e-^{-25). (3.7) 
v2 

Here M denotes the leading amplitude, including the CKM factor V^V^*, R e \,2 combinations 
of 5'C/(3)-breaking matrix elements, R u i,2 the dominant subleading contributions, and 5 the 
remaining one which is expected to be even smaller. The precise definitions of these parameters 
in terms of reduced matrix elements are given in the appendix. Finally we introduced for brevity 
A = A(l + A 2 /2), where A denotes the Wolfenstein parameter, and absorbed in these definitions 
the CKM factor R u = \/ p 2 + f] 2 and other numerical factors. Allowing, as discussed above, for 
SU (3) breaking of up to J"su(3) = 40%, see also appendix [Al and an enhancement of the penguin 
matrix elements such that they = 50% ofthe leading amplitude, we expect the 

order of magnitude of the parameters to be |i? e i.2|, |-R«i,2| ^ 10 — 20%, and \S\ < 1%. These 
values correspond to the combination of the generic scales rsu(3)> r pen with numerical factors like 
the CKM factor R u and Clebsch-Gordan coeffecients, given in the appendix. For 5, we assumed 
the dynamical suppression to amount to an additional factor of ~ 20% for the estimate. 

When discussing indirect CP violation, the B mixing phases <j> and <j> a enter. The latter 
is already constrained to be very small, and will be constrained more effectively by the corre- 
sponding golden mode B s — > J/ip(j), as it appears here only in suppressed decays. Therefore, 
we set this phase to its SM value, which affects the predictions for the indirect CP asymmetry 
in B s —7- J/^Ks only marginally. While ultimately these 13 (14) parameters should face 16 
measurements of decay rates and CP asymmetries, of these only 11 are available so far. The 
LHCb experiment is expected to determine the CP asymmetries in B s — > J/ifjKg soon, albeit 
with limited precision at first, yielding two more constraints. Therefore, for the time being, 
additional assumptions are necessary. It is however possible to make them in a transparent and 
testable way, which is one advantage of the present analysis. 

One obvious choice is the neglect of the A/m^-suppressed contributions discussed above. A 
qualitative confirmation of this assumption is given by decays dominated by the same topologies 
(but not related by symmetry), as e.g. B° — > D°(f>, to which an upper limit of ~ 10 -5 exists [21] , 
despite an enhancement factor of IV^/V^I 2 compared to that contribution in B° — > J/ijrir. 
Apart from this, here we also have the opportunity to test the contributions within the framework 
itself, which will be discussed in the next section. Making this assumption yields for the present 
situation 11 parameters for 11 measurements. 

Additionally, given that our main interest is the extraction of the mixing phase, we can use 
input for the CKM phase 7 to remove another free parameter. Unless otherwise stated, we 
will use 7 = (67.1 ± 4.3)°, extracted from a global fit [28 We will examine the possibility to 



3 Although S(B — > J/ipKs) enters this analysis, it does not affect the extraction of 7 significantly. We expect 
this input soon to be replaced by a competitive measurement in tree-level processes. 
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determine this angle within the fit below. These options show that the method proposed here 
is already applicable now and will be even more fruitful in the future. 

Considering non-SM contributions to this class of decays, we have two possible contribu- 
tions, which in general will both be present: new contributions to the mixing and to the decay 
amplitude. The former may be parametrized as (ft — > (ftsM + A^np, and result in an extracted 
value for the mixing phase (ft different from its SM value. This way to extract (ft takes into 
account the penguin "pollution", thereby excluding a source for a fake signal, and could then 



be used in an analysis including further observables from B mixing, as e.g. 129 30 . A model 



independent parametrization of NP in the decay amplitudes yields too many parameters to be 
fitted in full generality. While for a specific model flavour structure and relevant operators are 
usually given, allowing for the corresponding fit, in a model-independent approach we are forced 
to consider scenarios in which certain contributions are dominant, transforming in a definite way 
under SU(3). For minimal flavour violating scenarios |31| - |34| however, the parametrization given 
above remains valid, only the interpretation of the extracted parameter values changes, because 
the relative influence of Wilson coefficients and corresponding matrix elements is different. 



4 Phenomenology 

We now proceed with the analysis of the data shown in Tab. [TJ First we introduce combinations 



to 



of observables which obtain their main contributions from fewer parameters, similarly to 35 
understand better their physical interpretation. Then we start fitting the data, beginning with 
the SU(3) limit, to examine the necessity of the breaking terms introduced above. As there 
are less parameters involved here, we can allow for deviations from the limit 6 = 0, to check 
the compatibility of the data with our assumption. This is followed by an analysis including 
SU(3) breaking, where this contribution is again set to zero; here we first test the limit of 
vanishing penguin amplitudes, before allowing for both classes of contributions. We continue 
with possible future scenarios, to illustrate the impact of coming measurements, and finally 
discuss NP contributions to these decays to some extent. 



The software used for the fits is the "NLopt" software package 36 , together with the aug- 



mented lagrangian 37,38 and Sbplex/Subplex algorithms [36,39 . The best fit values have been 
checked additionally with Mathematica. 

In order to see what to expect from the fit, we introduce a common power-counting param- 
eter £ ~ 0.2 for the appearing quantities, where the assignment corresponds to rounding the 
logarithms log^(Ri,5), in accordance with the estimates given earlier in Sec. [3j 

AJ^U^I-^, and |<5|~£ 3 . (4.8) 

The combinations of observables facilitating the interpretation of the introduced paraemters are 
collected in Tab. [2] together with the leading order expressions using the power counting defined 
above. Estimates of parameters via these relations obviously only hold if the power counting is 
respected by the fit. 

The two isospin-related decay pairs offer the possibility to access the parameters R u i,2- the 
direct CP asymmetries are to leading order given by the corresponding imaginary parts, and 
the rate asymmetries Af, defined above, as well as 

2f (B° -> J/^vr ) - f(B- -» J/V.7T-) 
1 2T{B° J/Vw°) + T(B~ -> J/iftTT-) { ' ' 
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allow to constrain the real part of the combination AR U = R u \ — R u2 . As can be seen from 
Tab. [2j the former constrain the imaginary parts to be (at most) of the expected order, while 
the rate asymmetries yield a less consistent picure: Af prefers Re(Ai? u ) to be an order of 
magnitude larger than expected and negative (however with a significance of less than two 
standard deviations), while AJ in dataset 2 implies a very small value of at most the expected 
order, and dataset 1 yields a somewhat larger value, however with a positive sign. Combining 
these observations shows a slight inconsistency between the data and our theoretical expressions, 
specifically when considering the sum Af + A 2 AJ, which is expected to be safely below one per 
mill. Experimentally, its central value is several percent, dominated by Af, but enlarged in 
dataset 1 by the fact that AJ exhibits the "wrong" sign. If this is confirmed to higher precision 
in the future, it might become an actual hint of NP, making the corresponding measurements 
very interesting. 



Observable 


LO expression 


Experiment 


Af 


A 2 cos 7 [Re(i? M i) -Re(R u2 )} 


-0.037 ±0.025 


AJ 


- cos 7 [Re(R ul ) - Re(R u2 )] 


-0.13 ±0.06 






-0.01 ±0.05 


A C p(B° -> J/if>K°) 


2A 2 sin7lm(i? M i) 


0.010 ±0.012 


A CP (B- -> J/^K-) 


2 A 2 sin 7 Im( J R u2 ) 


0.001 ±0.007 


A C ?{B° -»• J/lpTT°) 


— 2sin7lm(i? u i) 


0.10 ±0.13 


A C p(B~ -> J/lpTT~) 


— 2sin7lm(i? U 2) 


0.01 ±0.07 






0.005 ±0.029 


AA CP 


4sin7lm(i? u ii? (: 2) 






— 2A 2 sin 7 cos(</>) Re(R u i) 




S(B -> J/^tt) + S(B -> J/ipK) 


2 sin 7 cos(<p) Re(R u i) 


-0.27 ±0.15 


R-kk — A 2 


-2A 2 [cos7Re( J R u2 ) + Re(3i? £ i + R e2 )\ 


-0.006 ±0.004 






-0.016 ±0.001 


Rkk - A 2 


-2A 2 [cos7Re(i?„i) + Re(3i? £ i - R e2 )] 


-0.015 ±0.004 



Table 2: Combinations of observables to disentangle the contributions of the different parameters, to- 
gether with their leading order (LO) expressions using the power counting discussed in the text. When 
two experimental values are present, the one for dataset 2 is given in the second line. 

The S'[/(3)-breaking parameters R e i, 2 are best determined by branching ratio combinations. 
The decay rate difference 

Kk = T(B S -+ J/^°) + 2T{B° ~ 2Re( ^ 2) (410) 

basically determines Re{R e2 ) to be at most of the expected order, while R^k = T(B~ — ► 
J/ifj-iT-)/T(B- -> J/i/jK-) and R KK = T(B S -»• J/i)K°)/F(B° -»• J/^K°) then mainly con- 
strain Re(R e i). An important point is that the value of Rkk (and the one of R^k in dataset 
2) can easily be explained by R € %, while the assumption of exact SU(3) symmetry would imply 
a very large real part of R u \ (R u2 for R^k), as their relative coefficient is 3/ cos 7 ~ 8. Re(R u \) 
determines at leading order in our expansion the differences 

AS X =Vx S(B° ^ X) +sm(j), (4.11) 

we are aiming to determine. Here it becomes therefore obvious that assumptions regarding SU (3) 
breaking influence the extraction of these differences strongly, which is one of the motivations 
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for the model-independent treatment of these terms in the present analysis. The mixing phase 
can be estimated from the relation 

- (VJ^KS(B d -> J/^Kg) + \ 2 Vj/i>7T S(B d -> J/^vr ) ~ (1 + A 2 ) sin0 = P 0.715 ± 0.023 , (4.12) 

where r\x denotes the CP eigenvalue of the final state X, and the main corrections cancel, as 
can be seen again from Tab. [2j 

A problem from the fitting point of view lies in the circumstance that none of the measured 
observables shown in this list depends on the imaginary parts of R e \^ on leading order. For 
that reason these imaginary parts are not well determined by the fit. For the same reason, 
however, predictions of observables will not be strongly influenced by them. In order for these 
parameters not to make the fit unstable, we constrain them to lie in a reasonable range, which 
we conservatively choose to correspond to rgu^\ = 40% of the leading matrix element, i.e. 
|Im(i? e i)| < 0.1 and |Im(J2 e 2)| < 0.14, which should affect only unphysical solutions. Another 
implication of this observation is, that these parameters do not effectively serve to fit the data, i.e. 
they should not be taken into account when determining the number of degrees of freedom in the 
fit. Here we simply reduce the number of parameters by two, which is again an approximation, 
but a reasonable one, given the situation that even huge values of these parameters are allowed. 
To remind the reader of this, we speak below of effective degrees of freedom. This situation could 
improve if the combination 

AA CP = A CP (B S -> J/iPK°) - A CP (B° ->• J/^vr°) (4.13) 

was measured to high precision, see once more Tab. [2j 

4.1 Fit in the SU(3) limit 

Performing the fit in the SU{3) limit (-R^i^ — * 0), and with 5 = yields a rather bad fit with 
Xmin/d.o.f. = 22.3(23.9)/5 for dataset 1(2), even when allowing up to l-R^i^l = 1- The reasons 
can be understood from the last paragraph, as the different estimates contradict each other in 
this limit. The situation worsens once reasonable values for the subleading matrix elements are 
enforced. Allowing for up to r pen = 50% yields x^/d.o.f . = 25.9(58.5)/5, where the latter 
result shows the strong influence of R^k- The outcome of such a fit is shown exemplarily in 
Fig.jlJ where the fit result for ASj/^k s versus S(B — > J/-07T ) is plotted; we emphasize that this 
is just for illustration. However, this plot allows for two connected interesting observations: the 
predicted sign of the penguin-induced phase shift is opposite to that in |5|, and Sqp(B — > J/ipn ) 
is predicted to be much smaller than the present central value. Both results are due to the 
precisely measured branching ratios entering the fit and will change again once SU(3) breaking 
is taken into account. 

We also repeat the fit, taking factorizable SU(3) breaking into account via the corresponding 
phase space and form factors, as this is advocated sometimes in the literature. This does not 
improve the situation, but yields Xmin = 28.0(23.9) for the two datasets. These decays are 
known anyway to be badly described by factorization, which is why we will not apply any such 
factors in the following. 

We proceed to investigate the influence of the parameter 5 neglected so far. Its inclusion 
does not improve the situation, either: the new fit yields a Xmi n /d.o.f. = 18.2(20.4)/3, which 
is worse than with (5 = 0. We take this observation as confirmation that the limit 5 — > is 
reasonable, and use this approximation in the following. 
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Figure 1: Fit result for ASj/^ Ks ver- 
sus S(B d -> J/^vr ) in the SU(3) limit 
for dataset 2, (x mm = 23.9), allowing for 
arbitrary values of the penguin matrix ele- 
ments. The inner (outer) area correspond 
to 68(95)% CL. Note the experimental 2a 
range for S(B — > J/tpn ), indicated by the 
light yellow (light grey) box. 



4.2 Fits including SU(3) breaking 

Including SU(3) breaking, the situation changes significantly. In a first step, we assume the 
subleading contributions ~ V* b Vcb to vanish identically. Obviously this implies vanishing direct 
CP asymmetries and the indirect ones to be equal to ± sin 0, such that the two measurements are 
effectively averaged to sin^ = 0.671 ±0.022. This fit works quite (very) well for the two datasets, 
we obtain x m in = 9.4(6.0) for 7 effective degrees of freedom. Not to count the imaginary parts 
of R e i as parameters is supported by the fact that setting them to zero does not change these 
values. We show the results in Fig. [2] therefore in the Re(i? e i) — Re(i? e 2)-plane, allowing the 
imaginary parts to reach values corresponding to 40% of the leading order matrix element. As 
can be seen there, both fits allow additionally for Re(i? e 2) = 0, while the best fit corresponds 
to Re(i? e i) = 0.047(0.059) respectively, yielding a value for the corresponding combination of 
S'C r (3)-breaking matrix elements of 19(24)% of the leading order one, which is perfectly within 
expectations. Therefore the data, strictly speaking, do not call for the inclusion of the penguin 
contributions at the moment (at least for dataset 2), although of course they are not perfectly 
described without them. In any case more precise measurements of the relevant CP asymmetries 
will decide the size of these terms in the future. 

Including again the parameters R u i,2, we obtain x m in = 2.8(2.3) for 3 effective degrees of 
freedom for the two datasets, when we refrain from applying strong restrictions on the parameter 
values^} For this fit, the parameters R € \ : 2 allow to accommodate the pattern of branching ratios, 
while the penguin contributions are mainly determined by the CP and isospin asymmetries. The 
central values of the absolute values l-R^i^l still tend to larger values than theoretically expected. 
This is not surprising, given the discussion on Aj ,7r above. The fit shows explicitly that the 
remaining Xmin ^ s mainly due to the branching ratios in B — > J/tpK: the fitted central values 
are one standard deviation higher (lower) for B° — > J/t()K (B~ — > J/ij)K~), underlining the 
importance of a new measurement of the ratio of these branching ratios, which correspondingly 
is predicted to take a significantly different central value than the one presently measured. 

Restricting the fit parameters to the expected ranges shows a difference between the two 

4 We do not allow for "exchanging roles" though, i.e. we continue to assume that TVo represents the leading 
matrix element. 
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Figure 2: Fit results (68/95% CL) in the Re(R el ) - Re(# e2 )-plane for datasets 1 (left) and 2 
(right), in the limit of vanishing penguin contributions. 



datasets: it does not worsen the fit as much as in the SU(3) case, but while for dataset 2 
it stays very close to the absolute minimum, Xminconstr = 2-8, it doubles for dataset 1 to 
Xmhfconstr =5.6. The new result for R w k obtained by LHCb seems therefore favoured by this 
fit. While it is too early to draw conclusions, this observation demonstrates once more the 
importance of precise branching ratio measurements in this context. The remaining difference 
stems from slight tensions in the B — > J/tpir CP asymmetries for dataset 1, which are both 
predicted to lie ~ la below their present central values. 

For both datasets, the shift ASj/^k now tends again to positive values, thereby lowering the 



corresponding tension in the UT fit 10-14 ; it is however still compatible with zero, in agreement 



with the observation above of a reasonable fit without penguin terms. The obtained ranges read 

AS j%k = [0.001, 0.005]([-0.004, 0.011]) , and (4.14) 
ASff^ K = [0.004, 0.011]([-0.003, 0.012]), (4.15) 

for 68% (95%) CL, respectively. The reason for the sign change compared to the SU(3) limit 
lies in the indirect CP asymmetry in B — > J/ip7r°, whose absolute value was forced to be smaller 
than the one in the golden mode before (see Fig. [TJ, but is now allowed to drive the shift to 
positive values. Differently put, neglecting SU(3) breaking drives ASj/^k to relatively large 
values with a different sign. S(Bd — > J/ipir°) is however still predicted by the fit to lie below the 
present central value of the measurement, thereby supporting the Belle result [40] over the BaBar 
one, which indicates a very large value for this observable [41]. These findings are illustrated in 
Fig. [3| The mixing phase is extracted as (j>^ t = 0.74 ± 0.03 (for both datasets), which is to be 
compared with 4*sM e = 0.73±0.03 when using the naive relation without penguin contributions. 
The inclusion of the correction therefore yields the same precision, but a slightly different central 
value. 

In Fig. [4] the prediction for the CP asymmetries in B s — > J /tpKs is shown for the two 
datasets. The allowed range remains rather large, which however indicates that the correspond- 
ing measurement will add important information to the fit already with moderate precision. In 
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Figure 3: Fit results for datasets 1 (left) and 2 (right), for ASj/^k versus Scp{B° — > J/ipir ), 
including all available data. The inner areas correspond to 68% CL and 95% CL with rgu^ = 
40% and '"pen — 50%. The outer one is shown for illustration purposes, only, and corresponds to 
95% CL when allowing for up to r su ^ = 60% and r pen = 75%. The light yellow area indicates 
the 2-<T range of the S(B° — > J/ipn ) average, the dashed line its central value. 
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A CP (B S -> J /iff K s ) A CP (B X JliffKs) 

Figure 4: Predictions for the two datasets for the CP asymmetries in B s — > J /ipKs, including 
all available data. The colour code is identical to Fig. [3} 

Fig. [5] we plot the correlation between ASj/^k s an d S(B S — > J/if;Ks), which allows to read off 
what impact the measurement of only this quantity would have on the phase shift. It will be 
possible to exclude parts of the allowed range for ASj/^k s with such a measurement; however, 
to see the full impact of future measurements, an analysis including all expected improvements 
is needed, which we will show in the next subsection. 
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Figure 5: Fit results for the two datasets for ASj/^k versus Scp(B s — > J/ipKs), including all 
available data. The colour code is identical to Fig. [3j 



Although in principle the fit is sensitive to the CKM angle 7 as well (see [2j[8]), the extracted 
range from the present data is huge. Actually this is not surprising, given the situation that 
neglecting all corresponding terms in the amplitudes leads to a rather good fit already. That 
situation may change, once a measurement of a significantly non- vanishing direct CP asymmetry 
is performed. 

4.3 Future scenarios 

As a main motivation for the present analysis are the expected precision measurements for 
the decays in question, we analyze their potential impact on the extraction of the weak phases 
described in the last subsections. For that purpose we define three scenarios, in which we assume 
improved measurements for some of the observables. Scenario 1 is corresponding to the situation 
in a few years, with ~ 5 fb _1 of LHCb data. Another few years later we have scenario 2, with the 
first ~ 5 ab _1 of SFF data. Finally scenario 3 corresponds to a situation in the farther future, 
with ~ 50 ab _1 of SFF data and 100 fb _1 of Super-LHCb data. As detailed studies are available 
only for a few of the observables in question, and a detailed analysis is beyond the scope of 
this work, we simply scale for the remaining observables the present statistical uncertainties to 
correspond to the integrated luminosities of the scenario in question. Details to that procedure 
and the resulting inputs can be found in appendix [B} 

To be able to compare with the results above, and at the same time to estimate the remaining 
dependence on the theoretical assumptions used in those fits, we add to the plots below another 
area, showing the allowed range at 95% CL when huge values are allowed for the parameters, i.e. 
r pen , r£j7(3) < 200%. Again, this is for illustration purposes only, to demonstrate the decreasing 
influence of the theoretical assumptions with additional data. 

Starting with scenario 1, the numerical improvement for the uncertainties of the mixing phase 
4> and the shift AS is small, as can be seen Tab. [3j However, Fig. [6] demonstrates the main effect 
of the new measurements: when removing the theoretical constraints, the range for AS does not 
increase much anymore, in contrast to the fit with present data, in which the upper limit grows 
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Figure 6: Fit results for scenario 1 for ASj/^k versus the mixing phase 4> (left), showing the 
three areas as in Fig. [3l and in addition the 95% CL range with r su ^ pen 



the right, we show the 95% CL regions of the three future scenarios, see Sec. 



and r 



pen 



50%. 



4.3 



200% allowed. On 
40% 



for 



r SU(3) 



by almost a factor of three. Although cancellations between some parameters remain possible, 
the physically relevant combinations are therefore under control with these measurements. 





Now 


Sc.l Sc.2 Sc.3 


5<t>(la) 

A(AS) (95%) CL 


0.030 
0.015 


0.027 0.021 0.015 
0.014 0.010 0.005 



Table 3: Results for the errors of the mixing phase <f> and the shift AS (A(AS') = AS max — AS'min) for 
the three future scenarios, comparing with the present situation. We do not show the central values, as 
they correspond to our choice of observables and are arbitrary in that sense. 

Regarding scenarios 2 and 3, the improvement in determining the mixing phase (ft stems from 
the one in S(B — > J/ipKs), and corresponds to the one that would be obtained using the naive 
relation for the mixing angle with AS = 0. 

Investigating the possible determination of the CKM phase 7, we observe that the allowed 
range remains very large (±35°), even within scenario 3. The reason for this is twofold: first 
of all the additional freedom introduced by the SU (3)-breaking parameters effectively prevents 
the penguin contributions to be determined precisely. Secondly, even with the precision for this 
scenario, some direct CP asymmetries are still compatible with zero, reducing the sensitivity 
to 7 further. This of course depends on the central values chosen, but the scenario is close to 
the present situation. However, to investigate the influence of this choice we repeated the fit 
with an idealized scenario, in which we have small SU(3) breaking and large imaginary parts of 
the penguin matrix elements. Even in that case we obtain only £7 ~ 15° for r pcn < 40%, and 
that range is only due to the fact that we chose parameters close to this limit, it is much larger 
without that assumption. 

We cannot confirm therefore the optimistic expected ranges for 7 given in p] within this 
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more general framework. Especially the estimate of S'C/(3)-breaking corrections given there 
does not seem general enough in light of the present analysis; the form factor ratio discussed in 
that context is not the main source of SU(3) breaking in these decays, as emphasized before. Of 
course, our theoretical understanding of SU (3) might improve during the next years; however, 
this kind of development seems hard to quantify. 



4.4 New physics contributions 

In the previous subsections we established working fits within the SM without a large enhance- 
ment for any hadronic parameter, therefore obviously no strong indication for NP in the decay 
amplitudes is observed. However, as an analysis within a Minimal Flavour Violation frame- 
work [3l] - |34| yields identical expressions, the possible enhancement of the subleading terms 
could as well be ascribed to NP contributions. In the latter case, however, additional contri- 
butions would be expected in processes with leading contributions from the same terms in the 
effective hamiltonian, such as B — > DD. 

Regarding NP in mixing, the extracted phase fits well with the value extracted from a UT 
fit with \V u b/V c b\ from semileptonic decays as main input for this quantity, yielding 0g M = 
0.84 ± 0.13, as does however the naively extracted value. 

The only actual tension is the one with B — > tv, which is measured well above the expected 



value [10 12- 14,22]. However, it is not clear if this tension should be discussed in the context of 
NP in mixing (note that the indication of NP in Bj mixing in [30] is mostly due to this mode, 
as there the assumption of NP residing only in mixing is made) : Firstly, this mode has a high 
sensitivity to NP, especially in form of charged scalars [42] . The very recent measurement of the 
branching ratios of B — > D(*'tv modes different from the SM predictions 43 strengthens this 
possibility. Secondly, for this tension to be a sign of NP in mixing (only) , the extracted values 
from leptonic and semileptonic decays should coincide, which they do only marginally. Finally, 
the experimental analysis for B — )• tu is extremely difficult. In any case, the observed shift due 
to penguin contributions is relatively small compared to that difference; it does therefore not 
offer an explanation for this tension. As mentioned before, it should however be included in the 
corresponding global analyses of NP in mixing. 



5 Conclusions 

The analysis presented here allows to include both, penguin and leading 5'C/(3)-breaking con- 
tributions to the golden mode, in a model-independent way. It thereby improves the extraction 
of the Bd mixing phase to match the impressive precision for S{Bd — > J/tpKg) of present and 
coming high luminosity colliders. The resulting uncertainty equals the one when using the naive 
relation with ASj/^ K = 0, see Fig. [6] 

The fits performed with present data show the importance of S'C/(3)-breaking corrections, 
without which neither a good fit nor a reliable value for the phase shift can be obtained. In 
fact, a reasonable description of the data can at the moment be achieved including only these 
corrections of O(20%), with vanishing penguin contributions. However, their inclusion yields an 
even better description of the data. The resulting phase shift is found to be small, \AS\ < 0.01, 
with a positive sign preferred, thereby slightly reducing the tension between different extractions 
of the CKM angle (3. This is the most precise determination of these contributions to date. 

An important observation is the importance of branching ratio measurements in this ana- 
lysis. Especially additional branching ratio measurements in B — > J /ipK (e.g. the ratio of the 
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charged and neutral modes) would be interesting, given the relatively large central value of the 
corresponding isospin asymmetry. Regarding the difference in the R n K measurements, the fit 
prefers the value obtained by LHCb, although this is not conclusive. 

The extraction of the CKM phase 7 turns out to be more difficult than expected. One 
reason lies in the fact that at the moment all direct CP asymmetries are compatible with zero, 
and, given the present data and projected experimental precision, several of them are likely to 
remain that way. However, even in an idealized scenario where they are larger, the precision in 
7 remains far from competitive without further knowledge of SU(3) breaking. Neglecting SU (3) 
breaking (or taking only the one by form factors into account), the pattern of branching ratios 
allows to determine the subleading terms much more precisely, which is why in that case a more 
optimistic result was achieved [8] . Of course, a better understanding of these effects might be 
achieved in the coming years, re-opening this possibility. 

The data do not indicate the presence of NP, although the penguin-induced shift in the 
mixing phase is too small to offer an explanation for the difference to B — >■ tv. However, the 
interpretation of this tension in terms of NP in mixing seems questionable, given its sensitivity 
to NP, and the situation in B — > D^*'tp and semileptonic decays. A general model-independent 
analysis of NP in the decay amplitude is impossible; for minimal flavour violating scenarios 
however the results shown here remain valid, and the enhancement of the penguin contributions 
compared to existing estimates might be attributed to NP contributions in that case. In this 
scenario, similar contributions should affect also e.g. B — > DD decays. 
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A Details of the SU(3) analysis 

We use conventions, in which the fundamental triplet is written as 

(u, d, s) = ((3) 1/2j i/ 2 ,i/3, (3)1/2,-1/2,1/3, (3)0,0,-2/3) • (A- 1 ) 

For the corresponding antitriplet the u has a negative sign, corresponding to the usual isospin 
representation. The meson states are then defined such that they correspond to "pure" SU (3) 
states without signs, i.e. every meson involving a u has an additional sign. This implies the 
assignments given in Tab. |4j 



M 


7T+ 


7T U 


7T 


m 


(IJz,Y) 


(1,1,0) 


(1,0,0) 


(1,-1,0) 


(0,0,0) 


M 


K+ 


K° 


K° 


K~ 


(I,Iz,Y) 


(1/2,1/2,1) 


(1/2,-1/2,1) 


(1/2,1/2,-1) 


(1/2,-1/2,-1) 



Table 4: Assignments of pseudoscalar octet meson states to SU(3) ones. 



The decompositions of the hamiltonian T-L u read, using the isoscalar factors given in 44 



lb ^ s 1,=, . /3. 



UV S ~ 2 (15)i + y-(15) + -(6)1 + y -(3) (A.2) 
15 



\(15)i/ 2 - J(6)i /2 + J|(3) 1/2 



(A.3) 



with AI 2 = 0, AY = —2/3 for the b — > s transition, and 

with AI 2 = —1/2, Ay = 1/3 for the b —> d one. Note that all occuring representations receive 
contributions from tree- and electroweak-penguin operators, while penguin operators contribute 
only to triplet matrix elements at this order. 



Decay 


(8|(3)|3) c 


B d -> J/^jK 


1 


B d -> J/ifjn 


1/V2 


B - -»■ J/V'-ftT" 


1 


B - -> J/^7T" 


1 


B s -> J/^vr° 





B s -»• J/V^K 


1 



<8|(3)|3) M 



(8|(6)|3) M (8|(15)|3) M 



1 

12 
i_ 

24 

1_ 

12 

1_ 

12 



Table 5: Coefficients for £? — > J/tpP amplitudes in the SU(3) limit. 

Projecting these operators on the relevant initial and final states yields the coefficients of the 
corresponding reduced matrix elements given in Tab.[5j Note that the amplitudes obey the trivial 
[/-spin relations for B~ — > J/ipP~ (P=ir,K) and between B° — > J/ipK° and B s — > J/ipK°, as 
well as isospin relations between the neutral and charged modes in B — > J /ipK and B — > J/ipir. 
We confirm the expressions obtained in 20 p 



Regarding SU (3) breaking, the tensor product of the leading-order hamiltonian T~L C with the 
isospin-conserving breaking term leads to 



H 



b— >s 



(8)0, 



0.0 



4( 15 )o,o, 



-2/3 + 2^ 3 -'o,o, 



-2/3 



and 



H 



1 



(8)o,0,0 ~ 4( 15 )l/2,-l/2,l/3 ~ y g( 6 )l/2,-l/2,l/3 - |( 3 )l/2,-l/2,l/3 ' ( A - 4 ) 

implying three additional reduced matrix elements, whose coefficients are given in Tab. [6| 

The corresponding coefficient matrix does not have maximal rank, i.e. it is possible to 
express the corresponding part of the amplitudes with less parameters, for example choosing the 
combinations A eX = (8|(3)|3)*- y/lJE <8|(15)|3>* and A e2 = (8| (6) |3>| + a/2/5 (8| (15) |3)^. Note 
that despite the somewhat small coefficient of the second matrix element in these relations, it 
is of course not clear which one will be dominating these combinations (if any). The dominance 
of one of the first ones would be signalled by only the corresponding one being sizable, while 
dominance of the absorbed one would result in — y2A e \ ~ A e2 in the following fits. 



For the parametrization given in Eq. (3.7) we introduced the following abbreviations: 
A/" = V c bV* s (8|(3)|3) c is the leading amplitude which is chosen to be real. The other matrix 
elements are normalized to the leading one, 

R el = J 4 el /(4(8|(3)|3) c ) and 



R e2 = A e2 /(2V2{8\(3)\3) c ) 



(A.5) 



5 In [20], a different normalization for the reduced matrix elements in A u is used. The translation reads 
|(15)|3)f = - > /i724(8|(15)|3> uI <8|(6)|3)f = <8|(6)|3) /2, and <8|(3)|3)f = ^3/8 <8|(3)|3\. 
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Decay 


(8|(3)|3)= (8|(6)|3): <8|(15)|3)« 


B d -> J/^K° 

B~ -> J/^K~ 
B- -> J/^vr- 


1 -i/i 

2 u Y 20 

FT i / i 

Y 32 4 Y 160 

1 o 

2 u y 20 

i fx FT 
4 Y 8 V 80 



i Fi FT 
4 Y 8 V 80 



Table 6: Coefficients for the linear S f C/(3)-breaking amplitudes from T~L C in B J/i^P decays. 



including the SU (3)-violating combinations of matrix elements introduced before, 



V33(8|(3)|3>„-(8|(15)|3>„ md v^8|(15)|3)„ 

v/l20(8|(3)|3} c *V15(8|(3)|3> C 

the dominant subleading contributions, and the remaining one 

^(8|(6)|3^-^(8|(15)|3) 
V60(8|(3)|3) c 



which is expected to be even smaller. Allowing for rgu{i) = |Ae1,2|/| (8|(3)|3) c | < 40% implies 
\Rei\ < 10% and \R e 2\ < 14%. Furthermore we impose in the constrained fit scenarios r pen = 



| (8|(N)|3) U |/| (8|(3)|3) c | < 50%, to be compared with r pen < 8% from |l6f|l8]. Note that this 
parametrization corresponds to the one given in [7] for the subclass of B u ^ decays, although 
the structure of the breaking terms looks different. The reason is that the £C/(3)-breaking term 
considered here, although conceptionally identical to one considered in that paper, includes 
a AU = part, which can be absorbed for [/-spin-related decays in the leading amplitude. 
Furthermore, the limit 5 — > corresponds to r\ = r 3 / 2 and = 3 / 2 in that paper. 



B Inputs for future scenarios 



We collect the estimated inputs for the different future scenarios discussed in Sec. 4.3 To find 
them, we proceed as follows: Explicit analyses are used where available; when they are absent, 
but a measurement is available from the corresponding experiment, the statistical uncertainties 
are scaled to correspond to the luminosity in the scenario considered, while leaving the sys- 
tematic errors untouched, i.e. cr now = J cr 2 at + cr 2 yst -> cr future = Jo^~L~Z /L future + of yst . 

We proceed the same way for £>-factory results to extract the expected SFF measurements, 
although this results in an even cruder estimate. Finally, when no measurement exists yet, we 
"guesstimate" the order of magnitude for the corresponding uncertainties, trying to stay on 
the conservative side. The listed uncertainties correspond to only one experiment, although of 
course in many cases different experiments will be able to measure them, thereby reducing the 
uncertainties further. Exceptions are the branching ratios in B — > J/ipK, for which we assume 
some progress despite them being already dominated by systematic errors, and the corrspond- 
ing CP asymmetries in scenario 1, where the number given corresponds to a guessed new world 
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Observable 



BR(B° - 
A CP (B° 
Scp(B - 
BR(B° - 
A CP (B° 
Scp(B - 
BR(B~ 
A C p{B- 

A CP (B- 

Rkk/% 
A C p{B s 
Scp(B s - 



> J/V^ u )/10" 4 
-> J/ipK°) 

> J/^7T°)/10- 4 
+ J/lpTT°) 

J/ipir°) 

> j/^k-)/io- 4 

-> J/'tpK') 

-> J/^7T-) 



Of, 



uture 



50Scl 5C?Sc2 5C?s. 



c3 



Ref. 



8.90 
0.012 
0.672 
0.178 
0.077 
-0.740 
9.68 
0.000 
3.78 
0.004 
3.67 
0.09 
0.01 



0.32 
0.021 
0.022 
0.016 
0.13 
0.15 
0.34 
0.007 
0.13 
0.029 
0.4 



0.25 
0.019 
0.020 
0.016 
0.13 
0.15 
0.25 
0.007 
0.08 
0.013 
0.25 

0.1 

0.1 



0.20 
0.015 
0.016 
0.008 
0.06 
0.06 
0.20 
0.003 
0.08 
0.012 
0.25 
0.1 
0.1 



0.20 
0.013 
0.012 
0.007 
0.021 
0.021 
0.20 
0.002 
0.07 
0.005 
0.22 
0.05 
0.05 



our guess (WA) 
our guess (WA) / [45] 
our guess (WA) / pB] 

our estimate 
chosen equal to Sep 



46 



our guess (WA) 
our estimate 
our estimate 
our estimate 
our estimate 
our guess 
our guess 



Table 7: Experimental data used for the analysis of the future scenarios 1-3. 0f utU re denotes the 
central value assumed for all scenarios. 

average. The results of this procedure are shown in Tab.[7j In order to create acceptable fits, we 
choose the projected central values to lie gaussian distributed around the values of an idealized 
scenario similar to dataset 2, with the uncertainties from future scenario 3. 
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